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Isolated neutron stars undergoing non-radial oscillations are expected to emit grav- 
itational waves in the kilohertz frequency range. To date, radio astronomers have 
located about 1,300 pulsars, and can estimate that there are about 2 x 10 8 neutron 
stars in the galaxy. Many of these are surely old and cold enough that their in- 
teriors will contain matter in the superfluid or superconducting state. In fact, the 
so-called glitch phenomenon in pulsars (a sudden spin-up of the pulsar's crust) is 
best described by assuming the presence of superfluid neutrons and superconducting 
protons in the inner crusts and cores of the pulsars. Recently there has been much 
progress on modelling the dynamics of superfluid neutron stars in both the Newto- 
nian and general relativistic regimes. We will discuss some of the main results of 
this recent work, perhaps the most important being that superfluidity should affect 
the gravitational waves from neutron stars (emitted, for instance, during a glitch) 
by modifying both the rotational properties of the background star and the modes 
of oscillation of the perturbed configuration. Finally, we present an analysis of the 
so-called zero-frequency subspace (i.e. the space of time-independent perturbations) 
and determine that it is spanned by two sets of polar (or spheroidal) and two sets of 
axial (or toroidal) degenerate perturbations for the general relativistic system. As in 
the Newtonian case, the polar perturbations are the g-modes which are missing from 
the pulsation spectrum of a non-rotating configuration, and the axial perturbations 
should lead to two sets of r-modes when the degeneracy of the frequencies is broken 
by having the background rotate. 

I. INTRODUCTION 



Jacob Bekenstein is one of those rare individuals who can make significant, original con- 
tributions to diverse areas of theoretical physics. He is also a man of great integrity and, I 
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believe, has a humility that serves him well in advising and supporting students and young 
scientists. I am profoundly grateful that fate allowed me to be one of those young scientists 
and now lets me participate in this celebration of his career. One of the areas of theoret- 
ical physics that Jacob has worked on is relativistic fluid dynamics. This is an important 
component of my current area of research, which is to develop models of Newtonian and 
general relativistic superffuid neutron stars. My original interest in superfluids, appropri- 
ately enough, was sparked by Jacob, when he suggested that I look at superffuid analogs of 
effects predicted for quantum fields in curved spacetimes (the Hawking and Fulling-Davies- 
Unruh effects). My current interest in superfluids is to determine how the dynamics of 
superfluid neutron stars differ from their ordinary, or perfect, fluid counterparts and if the 
different dynamics can lead to observable effects in gravitational waves. In the remainder 
of this article, I will give an overview of what my collaborators and I have accomplished 
so far, including some new results (from work with Nils Andersson) on the structure of the 
so-called zero-frequency subspace (i.e. the space of time-independent perturbations). The 
main purpose is to show that superfluidity in neutron stars should affect their gravitational 
waves in two ways, by modifying the rotational properties of the background star and the 
modes of oscillation of the perturbed configuration. 

While there are many mysteries about neutron stars that remain to be explained, we do 
have some significant observational facts to work with. For instance, Lorimer OJ reports that 
nearly 1300 pulsars (i.e. rotating neutron stars) have now been observed. By extrapolating 
the data on the local population, he can estimate that there are about 1.6 x 10 5 normal 
pulsars and around 4 x 10 4 millisecond pulsars in our galaxy. Of course, there are also 
neutron stars that are no longer active pulsars. To get a handle on their number Lorimer 
takes the observed supernova rate, which is about 1 per 60 years, and the age of the universe 
to find about 2 x 10 8 neutron stars in the galaxy. The overwhelming majority of these 
objects must be very cold in the sense that their (local) temperatures are much less than 
the (local) Fermi temperatures of the independent species of the matter. One can estimate 
the Fermi temperature to be about 10 12 K for neutrons at supra-nuclear densities, and it is 
generally accepted that within the first year (and probably much sooner than that) nascent 
neutron stars should cool to temperatures less than 10 9 K. This is an interesting fact, in 
that nuclear physics calculations of the transition temperature for neutrons and protons 
to become superfluid and superconducting, respectively, consistently yield a value that is 



3 



10 9 K in order of magnitude (for recent reviews see 0, |3|]). Thus we can expect that a 
significant portion of the neutron stars in our galaxy will have at least two (and perhaps 
more) superfluids in their cores. 

In addition to nuclear physics theory and experiment, the well-established glitch phe- 
nomenon in pulsars (e.g. Vela and Crab) is perhaps the best piece of evidence that 
supports the existence of superfluids in neutron stars. A glitch is a sudden spin-up of the 
observed rotation rate of a neutron star, and can have a relaxation time of weeks to months 
H . Baym et al have noted that a relaxation mechanism based on ordinary fluid viscosity 
would be much too short to explain a weeks to months timescale and so they argue that this 
signals the presence of a neutron superfluid. Now, a mainstay idea for explaining glitches is 
that of superfluids and their vortex dynamics, i.e. how the vortices get pinned, unpinned, 
and then repinned || ^|, [nj to nuclei in the inner crusts of the glitching pulsars. This is 
known as the vortex creep model and in it glitches are a transfer of momentum via vortices 
from one angular momentum carrying component of the star to another. The model has 
worked well to describe both the giant glitches in Vela and the smaller ones of the Crab. The 
vortex creep model can also be used to infer the internal temperature of a glitching pulsar, 
and for Vela it implies a temperature of 10 7 K |1(J. It is also interesting to note the work of 
Tsakadze and Tsakedze [1 f ] who have experimented with rotating superfluid Helium II and 
find behaviour very much like glitches in pulsars. 

The classic description of superconductivity in ordinary condensed matter systems is 
based on the so-called "BCS" mechanism (see, for instance, [|l2j or for excellent presen- 
tations): the particles that become superconducting must be fermions, and below a certain 
transition temperature there must be an (usually effective) attractive interaction between 
them (at the Fermi surface with zero total momentum). The interaction leads to so-called 
Cooper-pairing where a pair of fermions act like a single boson and a collection of them 
can behave as a condensate. The mechanism is very robust, which is why it also forms 
the basis for discussion of nucleon superfluidity and superconductivity Jl4j]; i.e. nucleons are 
fermions and the effective interaction between them at nuclear and supra-nuclear densities 
can be attractive. For instance, it is known experimentally that the lowest excited states 
in even-even nuclei are systematically higher than other nuclei because of pairing between 
nucleons which must be broken 



After many years of development, beginning with the work of Migdal flT5|] , a consistent 
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picture has emerged (in part, from gap calculations 0, 0]): At long-range the nuclear force 
is attractive and leads to neutron "Cooper" pairing in 1 So states in the inner crust, but 
because of short-range repulsion in the nuclear force and the spin-orbit interaction neutrons 
pair into 3 P2 states in the more dense regions of the core fIB |. In the crust protons are locked 
inside of neutron rich nuclei embedded in a degenerate normal fluid of electrons. In the inner 
crust the nuclei are also embedded in, and even penetrated by, the superfluid neutrons. In 
the core, however, the nuclei have dissolved and the protons remain dilute enough that they 
feel only the long-range attractive part of the nuclear force and pair in 1 So states. There is 
no pairing between neutrons and protons anywhere in the core since their respective Fermi 
energies are so different. The core superfluid neutrons and superconducting protons are also 
embedded in a highly degenerate normal fluid of electrons. Other possibilities, such as pion 
or hyperon condenstates, have been put forward but we will keep to the simplest scenario 
that considers only superfluid neutrons, superconducting protons, crust nuclei, and normal 
fluid electrons. 

There are several ways in which the dynamics of a superfluid differ from its ordinary 
fluid counterpart, and each difference should have some impact on the gravitational waves 
that a superfluid neutron star emits. One key difference is that a pure superfluid is locally 
irrotational. A superfluid, however, can mimic closely ordinary fluid rotation by forming a 
dense array of (quantized) vortices. In the core of each vortex the superfluidity is destroyed 
and the particles are in an ordinary fluid state, and can carry non-zero vorticity. A sec- 
ond, very important difference is when there are several species of matter in a superfluid, 
or superconducting, state. The superffuids of all the species will interpenetrate and each 
superfluid will be dynamically independent having its own unit four- vector and local particle 
number density. Lastly, superfluids have zero viscosity, but when vortices and excitations 
are present, then dissipative mechanisms can exist. For instance, the scattering of excita- 
tions off of the normal fluid in the vortex cores can lead to dissipative momentum exchange 
between the excitations and the superfluid, the net effect being that the superfluid motion 
becomes dissipative. This form of dissipative mechanism is known as mutual friction. 



In neutron stars there is a very efficient form of mutual friction [II, [IT], [L8|] that depends 
on the entrainment effect [|1^, 20], which Sauls H] describes as follows: even though the 
neutrons are superfluid and the protons are superconducting both will still feel the long- 
range attractive component of the nuclear force. In such a system of interacting fermions 
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the resulting excitations are quasiparticles. This means that the bare neutrons (or protons) 
are "dressed" by a polarization cloud of nucleons comprised of both neutrons and protons. 
Since both types of nucleon contribute to the cloud the momentum of the neutrons, say, is 
modified so that it is a linear combination of the neutron and proton particle number density 
currents. The same is true of the proton momentum. Thus when one of the nucleon fluids 
starts to flow it will, through entrainment, induce a momentum in the other fluid. Alpar et al 
| 17| have shown that the electrons track very closely the superconducting protons (because 
of electromagnetic attraction). Around each vortex is a flow of the superfluid neutrons. 
Because of entrainment, a portion of the protons, and thus electrons too, will be pulled 
along with the superfluid neutrons. The motion of the plasma leads to magnetic fields being 
attached to the vortices. The mutual friction in this case is the dissipative scattering of the 
normal fluid electrons off of the magnetic fields attached to the vortices. 



There has been much effort put forward to develop Newtonian [21, p2l 23, G4L E5| and 



general relativistic formalisms |26|, [27], |28|, [29|, |30], [31], ^1, [33|, |34| for describing superfluid 
neutron stars. In the simplest, but still physically interesting, formalism one has a system 
that consists of two interpenetrating fluids — the superfluid neutrons in the inner crust and 
core and the remaining charged constituents (i.e. crust nuclei, core protons, and crust and 
core electrons) that will be loosely referred to as "protons" — and the entrainment effect 
that acts between them. In principle the model can be expanded to have more than two 
interpenetrating fluids (see for instance). As well a given superfluid can be confined to 
a distinct region in the star [pq| . In this way the proton fluid, say, can be made to extend 
out farther than the superfluid neutrons. This is a first approximation at incorporating the 
fact that the superfluid neutrons do not extend all the way to the surface of the star. 

Our primary goal is to show that superfluidity will affect gravitational wave emission from 
neutron stars. We will see that a suitably advanced, but plausible, detector will have enough 
sensitivity at high frequency to see modes excited during a glitch. With such detections we 
will be able to place constraints, say, on the parameters that describe entrainment. But 
in addition to studying glitches, we also need to analyze further the recently discovered 
instability in the r- modes of neutron stars |5S| . The instability is driven by gravitational 
wave emission (the CFS mechanism |39|, fH], |41|) and the waves are potentially detectable by 
LIGO II P2] , H . The conventional wisdom early on stated that mutual friction would 
act against the instability in a superfluid neutron star and thus effectively suppress the 
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gravitational radiation. But Lindblom and Mendell |45| have found that mutual friction is 
largely ineffective at suppressing the r-mode instability. However, there are many questions 
about the spectrum of oscillation modes allowed by a rotating superffuid neutron star and 
the analysis of instabilities is very likely to be much richer than the ordinary fluid case. 
Thus, another goal here is to lay some groundwork for a future detailed study of the CFS 
mechanism in superfluid neutron stars. Specifically, we will demonstrate that the zero 
frequency subspace is spanned by two sets of polar (or spheroidal) and two sets of axial (or 
toroidal) degenerate perturbations for the general relativistic system. Like the Newtonian 
case ||6|], the polar perturbations are the g- modes which are missing from the pulsation 
spectrum of a non-rotating configuration, and the axial perturbations should lead to two 
sets of r-modes when the degeneracy of the frequencies is broken by having the background 
rotate. 

Below we will alternate between discussions based on Newtonian gravity and those us- 
ing general relativity. Accuracy demands that a fully relativistic formalism be employed, 
however there are some questions of principle for which a Newtonian formalism can suffice. 
For instance, in determining the number of different modes of oscillation that a superfluid 
neutron star can undergo it is much more tractable to use the Newtonian equations. But, 
ulitmately there is the need for a general relativistic formalism. Newtonian gravity does not 
include gravitational waves and so one needs a fully relativistic formalism to get an accurate 
damping time of a mode of oscillation due to gravitational wave emission. Also, there is 
the well-known problem that Newtonian models do not produce reliable values for the mass 
and radius, in that, for a given central density, the predicted mass and radius in Newtonian 
models may differ considerably from those of general relativity. This is a crucial point since 
the superfluid phase transition in a neutron star is sensitive to density, as are the parameters 
relevant to entrainment, and the oscillation frequencies can depend sensitively on mass and 
radius. 



II. GENERAL RELATIVISTIC AND NEWTONIAN SUPERFLUID 

FORMALISMS 

Based on the preceeding discussion, our formalism for modelling superfluid neutron stars 
must allow for two interpenetrating fluids (i.e. the neutrons and "protons") and the entrain- 
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ment effect. This must be the case whether we are working in the general relativistic regime 
or the Newtonian limit. A general relativistic superffuid formalism has been developed by 
Carter and Langlois |26], |27], [3(| [31], [32| and their collaborators p8| |29[ |33 



Here an 



action principle will be outlined that yields the equations of motion and the stress-energy 
tensor for this system. Although a variational principle also exists for the Newtonian regime 
\1E |, we will briefly indicate how to get the Newtonian equations by taking the appropri- 
ate limit of the general relativistic equations (see |^6| for full details). In the same spirit, 
we will use a "slow" velocity approximation to motivate an expansion that can be used to 



incorporate existing models of the entrainment effect |36 



A. The General Relativistic Formalism 



In this subsection the equations of motion and stress-energy tensor for a two-component 
general relativistic superfluid are obtained from an action principle. Specifically a so-called 



"pull-back" approach (for instance, see [^8|, |29[) is used to construct Lagrangian displace- 
ments of the number density four-currents that form the basis of the variations of the fluid 
variables in the action principle (they will also be used later in the analysis of the zero- 
frequency subspace). It will generalize to the superfluid case some of the techniques used to 
analyze the CFS mechanism in ordinary fluid neutron stars 0, |4l| . Finally, it can serve as 
a starting point for generalizing the Hamiltonian formalism developed by Comer and Lan- 
glois who limited their discussion to a superfluid in the Landau state [|l^] (i.e. purely 
irrotational). 



For both the rotation and mode calculations effects such as "transfusion" [33] of one 



component into the other (because of the weak interaction, for instance) will be ignored. 
The neutron and proton number density four- currents, to be denoted and p M respectively, 
are thus taken to be separately conserved, meaning 



V = 



(1) 



Such an approximation is reasonable for the mode calculations since the time-scale of the 
oscillations (which is milliseconds) is much less than the weak interaction time-scale in 



neutron stars, as long as the amplitudes of the oscillations remain small enough fl47| . In 
the case of slowly rotating neutron stars, it has been shown |48|, f49f that when the neutrons 
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and protons rotate rigidly at different rates then chemical equilibrium cannot exist between 
them. Of course, the energy associated with the relative rotation could be dissipated through 



a process like transfusion. But Haensel |50| has demonstrated that such a process would take 
months to years in mature neutron stars, and so again it will be neglected (since ultimately 
we are interested in gravitational waves emitted during a glitch, say, the timescale of which 
would be much shorter than the transfusion timescale). 
By introducing the duals to and p^, i.e. 

n uXr = e.ArX , n» = -^ vXt ti vXt , (2) 



and 



PuXt = e.ArX , P^ = jf vXT P^r , (3) 



respectively, then the conservation rules are equivalent to having the two three-forms be 
closed, i.e. 

V[^ Ar] = , V [a jw] = . (4) 

The reason for introducing the duals is that it is straightforward to construct particle number 
density three-forms that are automatically closed. The point is that the conservation of the 
particle number density currents should not — speaking from a strict field theoretic point of 
view — be a part of the system of equations of motion, rather they should be automatically 
satisfied when the "true" system of equations is imposed. 

This can be made to happen by introducing two three-dimensional abstract spaces which 
can be labelled by coordinates X A and Y A , respectively, where A, B, C etc = 1, 2, 3. By 
"pulling-back" each three-form onto its respective abstract space we can construct three- 
forms that are automatically closed on spacetime, i.e. let 

n vXr = N ABC (X D )V U X A V X X B V T X C , VvXr = P ABC (Y D )V U Y A V X Y B V T Y C (5) 

where N ABC an d Pabc are completely antisymmetric in their indices. Because the abstract 
space indices are three-dimensional and the closure condition involves four spacetime indices, 
and also that the X A and Y A are scalars on spacetime (and thus two covariant differentiations 
commute), the pull-back construction does indeed produce a closed three-form: 

V^rw] = V b (N ABC (X D )V U X A V X X B V T] X C ) = , (6) 
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and similarly for the protons. In terms of the scalar fields X A and Y A , we now have particle 
number density currents that are automatically conserved, and so another way of viewing 
the pull-back construction is that the fundamental fluid field variables are the spacetime 
functions X A and Y A [[Hj] . The variations of the three-forms can now be derived by varying 



them with respect to X A and Y A . 

Let us introduce two Lagrangian displacements on spacetime for the neutrons and pro- 
tons, to be denoted and respectively These are related to the variations 5X A and 
5Y A via a "push-forward" construction: 

5X A = -(y,X A )tt , 5Y A = -(y,Y A )^ . (7) 

Using the fact that 

VJX A = -v,([v M x A ]^) 

= - (v,x A ) v v e n - (v^x*) e n , (8) 



and similarly for the proton variation, then we find [33 



vXr 



SPuXt = ~ (^aPuXr + PaXr^ lAp + 2WV ' \£p + PuXa^T^p) = ~^ p PuXr , (9) 

where C is the Lie derivative. We can thus infer that 

W = P"V„# - e;V^ - p" (v ff £ + \g° p 5g^ . (10) 
By introducing the two decompositions 



n M = nu p , UfjU^ 1 = -1 , 



p» = pv» , VltV P = -1 , (11) 
then we can show furthermore that 

Sn = - V CT (tO - n (uvifVvil + \ W + u a u p ] 6<j 



crp J ) 
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Sp = -V, ( P C P ) - P (v u v°V a £l + \ [g ap + vV] 8g a ^ , (12) 



and 



<y„" = (5; + (^V„f£ - gV^') + ^'W^, . (13) 

We will take one more step ahead and associate a notion of Lagrangian variation with 
each Lagrangian displacement. These are defined to be 

A n = 5 + £^„ , A p = 5 + C ip , (14) 

so that it now follows that 

A n n pX r = , App^xr = , (15) 
which is entirely consistent with the pull-back construction. We also find that 

A n u» = X -u»u a u p A n g ap , A/ = W% p , (16) 



and 



A n dy\ Ta ^vXrag^^ ^ng fip > ^p^vXra ^vXrag^^ ' ' pp > (^^) 



A„n = ~ (g°o + u a u") A n g ap , A p p = -| (<7 CT " + v°v") A p g ap . (18) 



However, in contrast to the ordinary fluid case p0| , |41| , there are many more options to 
consider. For instance, we could also look at the Lagrangian variation of the neutron number 
density with respect to the proton flow, i.e. A p n, or the Lagrangian variation of the proton 
number density with respect to the neutron flow, i.e. A n p. It is not clear at this point 
how the existence of two preferred rest frames — one that is attached to the neutrons and 
the other that is attached to the protons — will affect an analysis of the CFS mechanism in 
superfluid neutron stars. 

Nevertheless, with a general variation of the conserved four-currents in hand, we can now 
use an action principle to derive the equations of motion and the stress-energy tensor. The 
central quantity is the so-called "master" function A, which is a function of all the different 
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scalars that can be formed from n M and p^, i.e. the three scalars n 2 = —n^n^, p 2 = —p^, 
and x 2 = —p^n^. In the limit where the two currents are parallel, i.e. the two fluids are 
comoving, then the master function is such that —A corresponds to the local thermodynamic 
energy density. In the action principle, the master function is the Lagrangian density for 
the two fluids. 

An unconstrained variation of A(n 2 ,p 2 ,x 2 ) with respect to the independent vectors n fl 
and p^ and the metric g^ v takes the form 

5A = fi^n" + XvW + I (™V + M , (19) 



where 



and 



fx,, = Brip + App , x M = Cp^ + An^ , (20) 



A dA B- 2 dA C- 2 dA (21) 

dx 2 dn 2 dp 2 

The momentum covectors /x M and Xn are dynamically, and thermodynamically, conjugate to 
n M and and their magnitudes are, respectively, the chemical potentials of the neutrons 
and the protons. The two momentum covectors also show the entrainment effect since it is 
seen explicitly that the momentum of one constituent carries along some mass current of 
the other constituent (for example, /i M is a linear combination of and p M ). We also see 
that entrainment vanishes if A is independent of x 2 (because then .4 = 0). 

If the variations of the four-currents were left unconstrained, the equations of motion 
for the fluid implied from the above variation of A would require, incorrectly, that the 
momentum covectors should vanish in all cases. This reflects the fact that the variations of 
the conserved four-currents must be constrained. In terms of the constrained Lagrangian 
displacements, a variation of A now yields 

5 (V=gA) = \v=g ( V + pV + nV) - ^V^g (W^C + p^x^Z) + t.d. 

(22) 

where the T.D. is a total divergence and thus does not contribute to the field equations 
or stress-energy tensor. At this point we can return to the view that n M and p^ are the 
fundamental variables for the fluids. Thus the equations of motion consist of the two original 
conservation conditions of Eq. ([!]) plus two Euler type equations 

^V^] = , p^V^Xu] = 0, (23) 
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since the two Lagrangian displacements are independent. We see that the stress-energy 
tensor is 

T? = *6Z + i?Xu + n' t ii v , (24) 
where the generalized pressure \I> is defined to be 



* = A - n'X - pt'xn ■ 



(25) 



When the complete set of field equations is satisfied then it is automatically true that 
V M 2? = 0. 

In a later section we will be interested in the linearized version of the combined Einstein 
and superfluid equations. It will thus be convenient to write down the variations of the 
momentum covectors in terms of the variations of the particle number density currents. 
Following the scheme of Carter pTfl , and Comer et al [|5"2"| , the variations of n P and Xp. due 
to a generic variation of and g pv take the form 



6fi p = A;6 Pa + B p °5n a + (S g A) Pp + (5 g B) 
5 Xp = C a 5 Pa + A a p 8n a + (5 g C) p p + (6 g A) ; 



n 



p ■ 



p • 



(26) 
(27) 



with 



A A -oHE _ 9 M _ 9 M _M 

P u •A-Qfiv ^^2^^P U Q^2^^^ u Flir3^^ v f)r r 2-P^^ l> ' 



B 



dp 2 
dB 



OA 



dp 2 
OA 



dx 2 



p" 



dC 



dn 2 



dA 



dx 2 



dA 



C pu - Cg pu 2 p p p v ^Q p2 P(p n iy) q x2 
and the terms 8„A, 6„B and 5„C are given by 



n p n v 



5 g A 



dA u „ dA u „ dA u v 
dn z dp z dx z 



&9. 



p v 



(28) 



(29) 



(5 g B and 5 g C being given by analogous formulas, with A replaced by B and C respectively). 



B. The Newtonian Limit 



The Newtonian superfluid equations can be obtained by writing the general relativistic 
field equations to order c°, where c is the speed of light, and then taking the limit that c 
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becomes infinite. To order c° the metric can be written as 

ds 2 = -c 2 (l + ^ dt 2 + bijd^dx 1 , (30) 

where the x % (i = 1,2,3) are Cartesian-like coordinates, and the gravitational potential $ 
is assumed to be small in the sense that — 1 << $/c 2 < 0. To the same order the unit 
four-velocity components defined earlier are given by 



u* = 1 - ^ + J- > «' = < , (31) 



c 



2c 2 ' 



and 

«* = l-*+||. « 4 = *J, (32) 



where v 2 = SijV^ v^ p , and and are the Newtonian three-velocities of the neutron and 
proton fluids, respectively. The three-velocities are assumed to be small with respect to the 
speed of light. The entrainment variable x 2 takes the limiting form 

,2 \ 



x 



2 



n n n p I 1 + ^ I , (33) 



where 

^ = (< - (< - $ (34) 
and we have introduced the new notation n n = n and n p = p (to distinguish between 
Newtonian quantities and their general relativity counterparts). Finally, we separate the 
"master" function A into its mass part and a much smaller internal energy part E, i.e. we 
write it as 

A = - [m n n n + m p n p ) c 2 - E{n 2 n , n 2 , x 2 ) , (35) 

where m n (m p ) is the neutron (proton) mass. 

To more closely agree with the Newtonian superfluid equations derived by other means 
p5| , a different choice for the independent variables is used which is the triplet of variables 
(n 2 , n 2 w 2 ). In this case E = E(n^, n 2 w 2 ) and the analog of the combined First and Second 
Law of Thermodynamics for the system takes the form 

dE = jj n dn n + n P dn p + adw 2 , (36) 

where 

dE dE dE 

on n on p aw 1 
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The generalized pressure to be renamed P, takes the limiting form 



P = -E + ji n n n + HpUp . 



(38) 



Formally letting the speed of light become infinite in the combined Einstein and general 
relativistic superfluid field equations, results in the following set of 9 equations: 







+ di (n n v l n j 
<)f +di(n p vi) 



Of 
dn, 



(39) 



and 



9 ( i 2a 



+vio j <+ 



2a 



^Si'O l ($ + — ) + 

m r . 



2a 
d 



5«5 kl (v l p - v l n ) d jV k n 
2a 



di 



+4d j [v; + 



2a 



VTLpUp 



+ 5 ij di $ + 



+ 



2a 



■5 ij 5 k i (v l n - v l p ) djv 



(40) 



The gravitational potential $ is obtained from 



did l § = AirG (m n n n + m p n p ) . 



(41) 



For this system having no entrainment means setting the coefficient a to zero. 

These equations are equivalent to those derived independently by Prix [p5|1 , using a New- 
tonian variational principle. We also note that they are formally equivalent to the two-fluid 
set developed by Landau |k| for superfluid He II. The two fluids in the Landau case are 
traditionally taken to be the normal fluid (i.e. the phonons, rotons, and other excitations) 
that carries the entropy and the rest of the fluid (i.e. the superfluid) that carries no entropy. 



C. An analytical equation of state with entrainment 

The item that connects the microphysics to the global structure and dynamics of super- 
fluid neutron stars is the master function, since it incorporates all of the information about 
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the local thermodynamic state of the matter. Ultimately, realistic models of superfluid neu- 
tron stars must be built using realistic master functions (i.e. equations of state). Only in this 
way can gravitational wave data be used to greatest effect to constrain the microphysics, 
such as parameters that are important for entrainment. Unfortunately, there is not yet 
a fully relativistic determination of entrainment (although Comer and Joynt are currently 
working on this), and the best that can be done is to adapt models that have been used 
in the Netwonian limit. Here we will describe an analytic expansion of the master function 



that facilitates the process (see [36] for all the details). 

In the limit where the fluid velocities are small with respect to c we see from Eq. ( |33"D 
that the combination x 2 — np is small. Thus, it makes sense to consider an expansion of the 
master function of the form [[36] 

oo 

A(n 2 ,p 2 ,x 2 ) = J2^(n 2 ,p 2 )(x 2 -np) 1 . (42) 

i=0 

The A, B, and C coefficients that appear in the definitions of the momentum covectors 
become 

i Xi(n 2 ,p 2 ) (x 2 — np" 

i=l 

ld\ 



A= — Xj(n 2 ,p 2 ) (x 2 — np] 

p l™dX if2 
-A > — [x -np 

n. n. r— : an \ 



B 



n dn n n r-f dn 
C = -^-HA--y^(x 2 -npY . (43) 

n f)r> r> n f)r> V / 



The expansion is especially useful for both the rotation and mode calculations, since x 2 = np 
for any zeroth-order, or background, quantity. In practice this means that only the first few 
Aj contribute, and for the mode calculations we need retain only Ao and Ai. The first 
coefficient Ao is directly related to equations of state that describe a mixture of neutrons 
and protons that are locally at rest with respect to each other. The other coefficient Ai 
contains the information concerning the entrainment effect. 
We will use a sum of two polytropes for A , i.e. 

X (n 2 ,p 2 ) = —m n n — a n n^ n — m p p — a v p^ v . (44) 

In Table | are given various parameter values that have been used in specific applications 



48, 52|. Model I corresponds to the values used in the initial study of modes on a non- 



rotating background by Comer et al |52], and is meant to describe only a neutron star core 
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without an outer envelope of ordinary fluid. Thus, the neutron and proton number densities 
vanish at the same radius. Model II is used in the mode study of Andersson et al |36| and 
is the most realistic. Its parameter values given in Table [I] have been chosen specifically to 
yield a canonical static and spherically symmetric background model having the following 
characteristics (see |36] for complete details): (i) A mass of about 1.4M , (ii) a total radius 
of about 10 km, (iii) an outer envelope of ordinary fluid of roughly one kilometer thickness, 
and (iv) a central proton fraction of about 10%. These are values that are considered to be 
representative of a typical neutron star. The background distribution of particles for this 
model is determined in Andersson et al |$6| and their graph and caption is reproduced here 
in Fig. |l|. Finally, model III is used by Andersson and Comer [48| to study slowly rotating 
configurations and is meant to be more realistic than model I by extending the radius out 
further and having a better total mass value, but with no distinct envelope so that the 
neutron and proton number densities vanish at the same radius. 

Our strategy for incorporating entrainment is to adapt models that have been used in 
Newtonian calculations. The particular model we will use is that of Lindblom and Mendell 
[ 15fl . It is a parametrized approximation of the more detailed model based on Fermi liq- 
uid theory ]55| developed by Borumand, Joynt, and Kluzniak ||56|| . As demonstrated in 
Andersson et al the Lindblom and Mendell model translates into the relation 

Ai = — ^— pnp , (45) 

Pnp PnnPpp 

where 

m n n = p nn + p np , m p p = p pp + p np , (46) 

and 

p np = -em n n . (47) 
Here e is taken to be a constant, which is the approximation introduced by Lindblom and 



Mendell []45 |. Prix, Comer and Andersson [^9[ argue that 



m p p 



m 



v 



m n n \m p 



-1 , (48) 



where m* is the proton effective mass. But, Sjoberg [[57]] has determined that 0.3 < m*/m p < 
0.8. Assuming a proton fraction of about 10% we can take as "physical" for a neutron star 
core those values that lie in the range 0.04 < e < 0.2. 



17 



III. SLOWLY ROTATING SUPERFLUID NEUTRON STARS 

The first application to consider is to the problem of axisymmetric, stationary, and asymp- 
totically flat configurations, i.e. the standard assumptions made for rotating neutron stars 
f58|, |59|, BO]. Although accurate codes exist for looking at rapidly rotating neutron stars in 



the ordinary fluid case [30, 61, 32 1, and can in principle be adapted to the superfluid case 
(Prix, Novak, and Comer, work in progress), we will limit our discussion to situations where 
rapid rotation accuracy is not needed. That is, we will describe models of rotating superfluid 



neutron stars that have been developed [|B, ^9], |63|] using a slow-rotation approximation. 



In the general relativistic regime, we discuss results of Andersson and Comer [IS| who have 



adapted the one-fluid formalism of Hartle |54| and Hartle and Thorne |55| to the superfluid 



case, whereas in the Newtonian limit it is the work Prix et al |49| (who have built on the 
Chandrasekhar-Milne approach |)6|, |67j) that will be reviewed. The most important aspect 



of the superfluid case in both the general relativistic regime and Newtonian limit is that the 
neutrons can rotate at a rate different from that of the protons, and this obviously has no 
analog in the ordinary fluid case. 

At the heart of the slow-rotation approximation is the assumption that the star is rotating 
slowly enough that the fractional changes in pressure, energy density, and gravitational field 



induced by the rotation are all relatively small p3 |. If M and R represent the mass and 
radius, respectively, of the non-rotating configuration, and Q n and Q p the respective constant 
angular speeds of the neutrons and protons, then "slow" can be defined to mean rotation 
rates that satisfy the inequality 

^or^or^ p «(|) 2 ^. (49) 

Since GM/Rc 2 < 1, the inequality also implies fl n , P R <C c, i.e. that the linear speed of the 
matter must be much less than the speed of light. 

This is not as restrictive as one might guess, especially for the astrophysical scenarios 
we have in mind. For a solar mass neutron star of radius 10 km, the square root of the 
combination on the right-hand-side of Eq. ( f49f ) works out to 11 500 s _1 . The Kepler limit 
(i.e. the rotation rate at which mass-shedding sets in at the equator) for the same star is 
roughly 7700 s^ 1 . The fastest known pulsar rotates with a period of 1.56 ms, which translates 
into a rotation rate of 4000 s _1 or nearly half the Kepler limit. This is why the slow-rotation 
approximation is still accurate to (say) 15 — 20 % for stars rotating at the Kepler limit. 
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This is seen in Fig. ^], which is taken from Prix et al [}49 |. It shows a one-fluid rotating 
star's equatorial (-R equ ) and polar (-R po ic) radii, determined using both the slow-rotation 
approximation and the very accurate LORENE code developed by the Meudon Numerical 



Relativity group |60| , |62|j . Notice that the slow-rotation approximation works quite well up 
to and including rotation rates equal to the fastest known pulsar. 

The slow rotation scheme is based upon an expansion in terms of the constant rotation 
rates of the fluids. The first rotationally induced effect that one encounters in general 
relativity is the linear order frame-dragging (i.e. local inertial frames close to and inside the 
star are rotating with respect to inertial frames at infinity). It is only at the second order 
that rotationally induced changes in the total mass, shape, and distribution of the matter 
of the star are produced. Since the frame-dragging is a purely general relativistic effect, the 
Newtonian slow rotation scheme yields nothing at linear order, but does have second order 
changes that parallel those of general relativity. 



Andersson and Comer j|8| apply their slow rotation scheme using the two-polytrope 
equation of state with the model III parameters listed in Table |. Although their formalism 
is general enough to allow for entrainment, they do not consider it in their numerical solutions 
to the field equations, since their main emphasis is to extract effects due to two rotation 
rates that can be independently specified. For instance, reproduced here in Fig. [5] are some 
of their typical results for the frame-dragging for modest values of the relative rotation 
Q n /Qp between the neutrons and protons. The solutions exhibit a characteristic monotonic 
decrease of the frame-dragging from the center to the surface of the star. Although likely 
very unrealistic, they also considered an extreme case of having the neutrons and protons 
counterrotate, with the result shown in Fig. (| (reproduced, again, from (^8|). Clearly the 
frame-dragging is no longer monotonic, and even changes sign. One can understand this 
behavior as follows: In the inner core of the star, the angular momentum in the protons 
dominates so that the frame-dragging is positive. But in the outer layers of the star the 
fact that the neutrons contain roughly 90% of the mass means they begin to dominate so 
that the frame-dragging reverses. Finally, we reproduce here one other result of Andersson 
and Comer in Fig. |^, which is the effect of relative rotation on the Kepler limit. When the 
relative rotation is larger than one there is little change in the Kepler limit, which is due 
simply to the fact that the neutrons contain most of the mass. On the other hand, when the 
relative rotation is being decreased toward zero, the Kepler limit rises because the frequency 
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of a particle orbiting at the equator is approaching the non-rotating limit (again because 
the neutrons carry most of the mass). 

Prix et al apply their Newtonian slow rotation formalism in a manner similar to 
Andersson and Comer. There are important differences, however, even beyond the exclusion 
of general relativistic effects, and these are (i) they use an equation of state that includes 
entrainment and terms related to symmetry energy (i.e. terms which tend to force the 
system to have as many neutrons as protons), and (ii) an exact solution to the slow rotation 
equations is used for the analysis. Thus, they are able to explore how entrainment and 
the symmetry energy affects the rotational configuration of the star. We reproduce here in 
Fig. |H their result for the Kepler limit as the relative rotation is varied, for different values of 
entrainment (denoted e) and symmetry energy (denoted as cr). The big surprise is that the 
symmetry energy has as much an impact as the entrainment. This fact has not been noticed 
before and should be explored in more depth, using a more realistic equation of state (like 

H). 

IV. THE LINEARIZED NON-RADIAL OSCILLATIONS 

The second application of our superfluid formalism is to the problem of non-radial oscil- 
lations. The ultimate goal is to calculate such oscillations, and the gravitational waves that 
result from them, for rotating neutron stars in the general relativistic regime. This is not 
an easy task, and the problem has not been solved fully (for rapidly rotating backgrounds) 
even for the "simpler" case of the ordinary perfect fluid. That is, while some recent progress 



has been made to calculate the frequency of oscillations |69| , there are as yet no complete 
determinations of the damping rates of the modes because of gravitational wave emission. 
Even using the slow rotation approximation there are complications, due to questions about 
the basic nature of the modes and if they can be separated into purely polar and axial parts, 
or if they are of the inertial hybrid mode class (as in |70| , [H|1). Nevertheless, we can gain 
valuable insight by considering non-radial oscillations on non-rotating backgrounds. We will 
use the Newtonian equations to reveal the nature of the various modes of oscillation, by 
given the highlights of the recent analysis of Andersson and Comer fl4q| , and we summarize 
the main results of general relativistic calculations [^, [52| of mode frequencies and damping 
rates. 



20 



A. Linearized oscillations in Newtonian theory 

The investigation of the nature of the modes of oscillation in both the Newtonian and 



is 



general relativistic regimes has now a decade and a half of history. Epstein's work [47 
the beginning, since he is the first to suggest that there should be new oscillation modes 
because superfluidity allows the neutrons to move independently of the protons, and thereby 
increases the fluid degrees of freedom. Mendell |22] reaches the same conclusion and moreover 
argues, using an analogy with coupled pendulums, that the new modes should have the 
characteristic feature of a counter-motion between the neutrons and protons, i.e. in the 
radial direction as the neutrons are moving out, say, the protons will be moving in, which is 
to be contrasted with the ordinary fluid modes that have the neutrons and protons moving in 
more or less "lock-step" . This basic picture has been confirmed by analytical and numerical 



studies |23|, |36|, [46], 0, [7^, [73|, |74J and the new modes of oscillation are known as superfluid 
modes. As we will see below they are predominately acoustic in nature, and have a sensitive 
dependence on entrainment parameters. It is worthwhile to mention again that our equations 
are formally equal to the two-fluid set of Landau for superfluid He II. With some hindsight, 
one recognizes that the superfluid modes could have perhaps been inferred to exist from 



a thermomechanical effect [TI| in which the normal fluid and superfluid are forced into 



counter-motion by passing an alternating current through a resistor placed in the container 
that holds the fluid. 

The existence of the superfluid modes seems to confirm one's intuition that a doubling of 
the fluid degrees of freedom should lead to a doubling of modes. A review of the spectrum of 
modes for the ordinary fluid should then give one an idea of what to expect in the superfluid. 
McDermott et al |75| have given an excellent discussion of many of the possible modes in 



neutron stars, and these include the polar (or spheroidal) f-, p-, and g-modes and the axial 



(or toroidal) r-modes. But a puzzling aspect in all of this is that Lee's |72| numerical analysis 
does not reveal a new set of g-modes, and in fact does not find any g-modes of non-zero 
frequency. Of course, the model he considers is that of a zero-temperature neutron star, 
and so one would not really expect g-modes like those of the sun (which exist because of 
an entropy gradient) to be important in a mature neutron star. However, Reisenegger and 
Goldreich |76| have shown conclusively that a composition gradient, such as the proton 
fraction in a neutron star, will also lead to g-modes, and the model of Lee does have a 
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composition gradient. Fortunately, this issue has been clarified by Andersson and Comer f46 
who use a local analysis of the Newtonian equations to confirm Lee's numerical result that 
there are no g-modes of non-zero frequency. Their analysis of the zero-frequency subspace 
reveals two sets of degenerate spheroidal modes, which they take to be the missing g-modes. 
They also find two sets of degenerate toroidal modes, which they interpret to be r-modes. 
And in fact when they add in rotation they find that the degeneracy is lifted and two sets 
of non-zero frequency r-modes exist. 



Apart from questions about the g-modes, Andersson and Comer |^6[ also illuminate the 
character of the superfluid modes. Although they do not solve the linearized equations for 
global mode frequencies, they are able to demonstrate that the equation that describes the 
radial behaviour of the superfluid modes is of the Sturm-Liouville form for large frequencies 
uj. Thus, one can expect there to be a set of modes for which — > oo as the index n — > oo. 
The equation that describes the ordinary fluid modes is also of the Sturm-Liouville form for 
large frequencies and so it also has a set of modes with the same mode frequency behavior. 
That is both sets of modes will be interlaced in the pulsation spectrum of the neutron star. 



Finally, Andersson and Comer [16 also use a local analysis of the Newtonian equations 



to find a (local) dispersion relation for the mode frequencies. Letting 



^2 ^"n 9fi n ^2 Tip 9 Hp (50) 

n m n dn n ' p m p dn p ' 

and assuming that the proton fraction (i.e. the ratio of the proton number density over the 
total number density) is small, then they find that one solution to the dispersion relation is 

ul « Ll (51) 

where 



, K '(' +iM . (52 ) 



Here c n is essentially the speed of sound in the neutron fluid, I is the index of the associated 
spherical harmonic Y™ of the mode, and r is the radial distance from the center of the star, 



and so this is the classic ordinary fluid solution in terms of the Lamb frequency L\ [77 
Likewise, they find another solution of the form 
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where c p is roughly the speed of sound in the proton fluid. Thus both solutions are of 
predominately acoustic nature, but we see that the second solution, which corresponds to 
the superfluid mode, has a sensitive dependence on entrainment (through the appearance of 
the proton effective mass). An observational determination of the superfluid mode frequency 
via gravitational waves, say, could be used to constrain the proton effective mass [^, |78f . 
This would translate into a deeper understanding of superfluidity at supra-nuclear densities 
since the effective proton mass is part of the input information for BCS gap calculations 
7|. 



B. Quasinormal modes in general relativity 

We have just seen that the spectrum of mode pulsations in a superfluid neutron star is 
significantly different from its ordinary fluid counterpart. Moreover, we have also seen that 
the superfluid modes have a sensitive dependence on entrainment. We will now confirm, and 
build on, this basic picture in a quantitative way by looking at numerical results for the modes 
obtained using the general relativistic formalism. The key results to be described come from 
the work of Comer et al and Andersson et al ||36|| . It is worth noting that the set of 



linearized equations that describe the mode oscillations of superfluid neutron stars has much 
in common with the ordinary fluid set, and so many of the computational and numerical 
techniques that have been developed for the ordinary fluid pOL |82L B3, |84L [|6L ET], |88fl 



can be adapted to the superfluid case. That being said, Andersson et al [|36] have developed 
a new technique for calculating the damping rates of the modes due to gravitational wave 
emission. 

Before we get to the ordinary fluid and superfluid modes, there is another set of modes, 
called w-modes, that we will discuss that exist only in a general relativistic setting. They 
were first discovered by Kokkotas and Schutz |89| and are due mostly to oscillations of 
spacetime, coupling only very weakly to the matter. For instance, Andersson et al pC§ use 



an Inverse Cowling Approximation where all the fluid degrees of freedom are frozen out and 



were able to find the w-modes. Comer et al |52| have obtained w-modes in the superfluid 



neutron star case and find that they look very much like those of ordinary fluid neutron 
stars. As well they do not find a second set of w-modes because of superfluidity. Both 
results are due to the fact that the w-modes are primarily oscillations of spacetime. 



23 



From the previous subsection we expect that there will be no g-modes in the pulsation 
spectrum (cf. Sec. |V| where we find two sets of polar perturbations in the zero frequency 
subspace), but there should be interlaced in it the ordinary and superfluid modes. In general 
relativity one obtains quasinormal modes because each frequency will have an imaginary 
part due to dissipation via gravitational wave emission. Such modes correspond to those 
particular solutions that have no incoming gravitational waves at infinity. In Fig. [7|, taken 
from []3~6|| , is graphed the asymptotic amplitude of the incoming wave versus the real part 
of the mode frequency for model I of Table |. The zeroes of the asymptotic incoming wave 
amplitude correspond with the deep minima of the figure. As expected, the ordinary and 
superfluid modes are interlaced in the spectrum and the lowest few have been identified in 
the figure. 

We recall that the ordinary fluid modes should be characterized by the neutrons and 
protons flowing in "lock-step" whereas the superfluid modes should have the particles in 
counter-motion. As well, the superfluid modes should have a sensitive dependence on en- 
trainment. Both are confirmed by the next two figures (Figs. |8] and ^|both taken from [p6fl), 
which also reveal a phenomenon known as avoided crossings. Fig. ^| graphs the real part 
of the first few ordinary and superfluid mode frequencies as a function of the entrainment 
parameter e (for the "physical" range discussed earlier in Sec. II). The solid lines in the 
figure are for the ordinary fluid modes, and we see that the first few are essentially flat as 
the entrainment parameter is varied, but the superfluid modes (the dashed lines) are clearly 
dependent on the entrainment parameter. Fig. |^ is a graph of the Lagrangian variations in 
the neutron and proton number densities. The two plots on the far left of the figure are 
for the modes whose frequencies are labelled by ct (for the superfluid mode) and b (for the 
ordinary fluid mode) in Fig. || The ao graph in Fig. || clearly indicates a counter-motion of 
the neutrons with respect to the protons, whereas the &o graph shows the opposite behavior. 

Another obvious feature of both figures is the avoided crossings phenomenon. Near the 
top of Fig. H we see that there are points in the (e, Re uM) plane where the solid and 
dashed lines approach each other, but just before crossing they diverge away from each 
other. The most interesting aspect of an avoided crossing is how the mode functions behave 
before, during, and after the avoided crossing. This is shown in Fig. ^, where the middle 
and right-hand-side graphs are for the modes of Fig. [8] labelled by {a 01 , 6 .i} and {a .2, &0.2}, 
respectively. In the middle graphs we see that the modes no longer have a clear distinction as 
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to whether or not the neutrons and protons are flowing together or in counter-motion. But 
in the {a .2, &0.2} graphs of Fig. || we see that now it is the superfluid modes that have the 
particles flowing together and the ordinary fluid modes show the counter-motion. Although 
we will not go into details here, Andersson and Comer and Andersson et al have 



suggested that this may explain one of the interesting results of Lindblom and Mendell ^5 
on the effect of mutual friction damping on the r-modes in superfluid neutron stars, and 
that is that mutual friction damping is negligible for the r-modes except for a small subset 
of the values of the entrainment parameter e. Mutual friction should be most effective when 
the neutrons and protons are in counter-motion, as in the superfluid modes, and what may 
be happening is that as the entrainment parameter is varied it is possible that they find 
modes beyond an avoided crossing where the ordinary fluid modes take on the characteristic 
of counter-motion. 



C. Detectable Gravitational Wave Signals? 

While the study of modes in superfluid neutron stars is a fascinating and complex mathe- 
matical and theoretical physics problem, the real hope is that one can use the results to make 
scientific progress. This is why it is important to understand superfluid neutron star dynam- 
ics for realistic astrophysical scenarios, because one wants to know if there are detectable 
gravitational waves that carry imprints of superfluidity. In other words, one wants to de- 
velop a gravitational wave asteroseismology as a probe of neutron star interiors, in much the 
same way that helioseismology is already a probe of the sun and asteroseismology is a probe 
of distant stars |)T|. This exciting possibility is being discussed |j78|, [93|, [95| in the 



literature and already some quantitative statements are in hand. Unfortunately, estimates 
for LIGO II suggest that one needs modes of unrealistic amplitudes to ensure detection. The 
best possibility is for detection of modes following neutron star formation from gravitational 
collapse, but even this has to be qualified |52, Q because of low event rates and uncertain- 



ties in the energy that will get deposited in the oscillations. However, there is no reason 
why one should only be pessimistic, since clearly gravitational wave detection will improve 
as more experience is gained, and new technology will also lead to improved sensitivities. 
For instance, there is already the so-called EURO [^6] detector being discussed, which is 
a configuration of several narrow-banded (cryogenic) detectors operating as a "xylophone" 
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that should allow high sensitivity at high frequencies. 



Andersson and Comer ]7S| and Andersson et al |36] have taken the estimated spectral 
noise density for such a configuration and have used it to determine the signal-to-noise for a 
detection of oscillation modes from superfluid neutron stars that have been excited during 
a glitch. These estimates have been for two of the best studied glitching pulsars, which are 
the Crab and Vela pulsars. One assumes that a typical gravitational wave signal from a 
mode takes the form of a damped sinusoidal, where the damping time is that of the mode 



itself |36|, |78| . The amplitude of the signal can be expressed in terms of the total energy 
radiated through the mode. Andersson and Comer and Andersson et al assume that this 
total energy is comparable to the amount of energy released during a glitch, which for the 
Crab and Vela pulsars can be of order 10~ 12 — lCU 13 M Q c 2 |97| , |98]. Unfortunately, an error 
in the signal-to-noise calculation of Andersson and Comer, to be corrected in Andersson 
et al, has incorrectly estimated the predicted signal-to-noise for the EURO configuration. 
Fortunately, the revised data still indicate sufficient sensitivity to expect detection for Crab- 
and Vela-like glitches. 



V. THE GENERAL RELATIVISTIC ZERO-FREQUENCY SUBSPACE 

Andersson and Comer j|6| have determined that the zero-frequency subspace in Newto- 
nian theory is spanned by two sets of polar and two sets of axial degenerate perturbations. 
These solutions are time-independent convective currents. They were stated to be the two 
missing sets of g-modes and the two sets of r-modes that become non-degenerate when ro- 
tation is added. We will extend this analysis to the general relativistic case and show that 
there are two sets of polar perturbations, thus supporting our earlier claim that there will be 
no non-zero frequency g-modes in the pulsation spectrum of a general relativistic superfluid 
neutron star. We will also find two sets of axial perturbations, which could presumably lead 
to two sets of r-modes (or more general hybrid modes) when the background rotates. 
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A. The background spacetime and fluid configuration 



The background is treated exactly as in Comer et al. ||52|| , i.e. it is spherically symmetric 
and static, and thus the metric can be written in the Schwarzschild form 



ds 2 = _ e »(r) dt 2 + e X(r) dr 2 + f 2 + _ 



(54) 



The two conserved currents n v and p u are parallel with the timelike Killing vector t v = 
(1, 0, 0, 0), and are thus of the form 

n u = n(r)U u , p v = vij)U u , (55) 

where W = t u /\t\. Likewise, the chemical potential covectors fi u and Xv become 

Xu = x(r)U u , fi v = fi{r)U u , (56) 



where fi = Bn + Ap and x = Cp + -An. Explicit solutions for the background configurations 
can be constructed following the procedure of Comer et al 



B. The linearized fluid and metric variables 



Making no assumptions yet on the metric and matter variations, we will first insert the 
background metric and matter variables into Eqs. (|I"2"D, (|1~3|) and fl2T|). Without too much 
effort, one finds 



for the fluid velocity perturbations, 
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for the density variations and 



p(lX + ^t + ™<) . (58) 



Spio = ^e^ 2 5goo-e^ 2 (B 5n + A° 5p 
5fa = ^ 2 6g 0t + e^ 2 9ij {Bn^i + Apjft 
5Xo = lxe- u/2 Sg 00 ~ e^ 2 (C° 5p + A° 5n) , 



S X i = Xe- u/2 Sg 0i + e-^g^ ( Cp J^j + An^{ ) , (59) 



for the momentum covector variations, where Aq, B®, and C® can be calculated from Eq. (f28[). 
One can also show that 

5 A = -fi8n - x Sp , 6V = (B%n + A° p) 5n + (c$p + A° nj dp . (60) 

Since the four velocity perturbations must be time-independent we see that the Lagrangian 
displacements must take the form 

£ = e^tfti 4 + £ , Q = e^HSv 1 + Q , (61) 

where and Q are "integration constants" (i.e. they are independent of time but depend 
on the spatial coordinates). We will see below that their only role is that once they have 
been determined they specify Sn and Sp via Eq. ( j58|) . 

Although the Einstein equations must be analyzed using a decomposition in terms of 
spherical harmonics for the perturbations, we can actually completely solve the linearized 
Euler equations, which take the form 

dtSfM = di5fj, t , d t 5xi = di5xt ■ (62) 

One can easily verify that the left-hand-sides of each equation is zero (by taking a time 
derivative of 8 fa and 5\i using Eq. (|59|) ) and it thus follows that 



6m = , 5 X t = 0. (63) 

From Eq. (|59| ) we see that these last two equations can be used to find 5n and bp in terms of 
5<7oo- Given that the conservation equations are satisfied automatically by the Lagrangian 
displacements, then all of the fluid equations have been solved. 
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We have seen earlier that the question of chemical equilibrium is an important aspect of 
the types of perturbations that can be induced on a superffuid neutron star. Langlois et al 
33| have argued that the general condition for chemical equilibrium to exist between the 



two fluids is that 

(5 = v v 0, - X u) = . (64) 

For perturbations that do not maintain chemical equilbrium then it is the case that 5/3 ^ 0. 
Using the variations above we find 

5(3 = (A° - £ °) 5n - (A* - C °) 5p . (65) 

We note here that axial perturbations are such that 5n = 5p = and thus it follows that 
5/3 = 0, that is axial perturbations on spherically symmetric and static backgrounds must 
necessarily maintain chemical equilibrium between the two fluids if the background fluids 
were in chemical equilibrium. 



C. The Zero- Frequency Subspace 



We have exhausted the information that can be extracted from just using the background 
configuration in the formulas for the variations. To finish mapping out the zero-frequency 
subspace we must examine the Einstein equations. In doing this it will be very convenient 
to consider two types of perturbations and those are the polar (or spheroidal) and the axial 
(or toroidal) perturbations. For the metric perturbations we will use the Regge- Wheeler 
gauge |8"0fl, in which the polar components of the metric perturbations can be written as 



e u ^H (r) 






fli(r) 
e x ^H 2 (r) 
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(66) 



and the axial components as 
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(67) 
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where the Y™ are the spherical harmonics. We write the polar unit four-velocity perturba- 
tions as 



-z//2 



I H 
±VJr)4- n 



and the axial perturbations as 



e -u/2 



-v/2 



6 A u 



r 2 sin# 





UJr)-g 
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(68) 



U P {r)ie 
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(69) 



Although we have already seen that the particle number density perturbations can be 
solved in terms of <5goo> it is worthwhile to comment just a little more on their form. The 
easy case is that of axial perturbations since for them 5 A goo = and so 5 A n and S A p must 
both vanish for a generic master function. The polar perturbations in the particle number 
densities are a bit more complicated to determine, but in the end take a simple form. Since 
5 p n and 5 p p must both be time independent, a setting of time derivatives of Eq. ([58]) to 
zero yields constraints on the velocity perturbation functions, i.e. 



l{l + l)V n 



11 



(rne x/2 W n )' , l{l + l)V p 



P 



rpe x/2 W p )' . 



(70) 



But recall that the Lagrangian displacements still have the "integration constant" terms. 
However, using the same type of decomposition for as used for 5 p u l and S p v 1 above, and 
in place of the W ntP (r) and V n)P {r) coefficients we substitute some new coefficients A ntP (r) 
and Bn tP (r), say, then it follows that 



5 p n = 5n(r)Y{ 



5 p p = 5p(r)Y{ 



(71) 



where the Sn(r) and Sp(r) are linear combinations of A np (r) and B n ^ p (r) and their deriva- 
tives. These new coefficients A n)P (r) and B n>p (r) appear nowhere else but in 5 p n and S p p 
and thus this is why we stated earlier that the only role of the "integration constants" is to 
determine the particle number density perturbations. Of course at this point we can forgo 
using A n<p {r) and B ntP (r) and just solve for the Sn(r) and 8p(r) instead. 
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After some algebra it can be shown that the perturbations for / > 2 result in three distinct 
groups of the linearized Einstein and superfluid field equations, and these are 



(i) I > 1 Group I: 



= e -\- 2 K" + e- x 3 \rK' 



'1(1 + 1) 
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|if = ^n + C >. 
(ii) Z > 2 Group II: 

167rre A ... ... . 

1 + z(T+T) (/U + xp p) ' 
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= e -(^ A )/ 2 (e^-^H,)' + 16ne x ([inV n + X pV P ) . (73) 
(iii) I > 2 Group III: 

= ^ — K + ( ^ e A J h - 167re A { XP U P + finU n ) , 

= (I- l)(l + 2)h 1 , 

= e-^' 2 (e^-Wh)' . (74) 

A quick counting of the number of independent functions, and the number of equations, 
shows that Group I appears to have more equations than unknowns. However, such a 
result is not unexpected because of the Bianchi identities. For the present discussion, the 
important point about the Group I equations is that their solutions represent a subset that 
map static and spherically symmetric stars, with no mass currents, into other (nearby) static 
and spherically symmetric stars, also having no mass currents. More interesting counting 
comes from the Group II and III equations. For Group II one can show that the last equation 
in the group is a consequence of the other three. Thus, we can specify arbitrarily W n and 
W p , for instance, and then the other variables (Hi, V n , and V p ) can be determined from 
the field equations. Likewise, for Group III we can specify freely U n and U p , and then the 
remaining variable ho is determined from its field equation (because it is clear that hi — 0). 

For various reasons, the cases of / = 0, 1 must be distinguished from that of I > 2. One 
reason is that there is more gauge freedom, which allows us to set K(r) = hi(r) = for 
/ = 1 and in addition Hi(r) = for I = (which also has no axial perturbations). Without 
listing all the formulas, we find that the counting of the number of equations and unknown 
functions is similar to I > 2. In particular, the I = 1 analysis of the Groups II and III 
equations reveals that each have two arbitrary functions that must be specified before a 
solution can be obtained. For 1 = 0, the polar currents must also vanish (otherwise they 
would diverge at the center of the star). Hence, there are only I = solutions that map 
static and spherically symmetric stars into other static and spherically symmetric stars. 
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D. Decomposition of the zero-frequency subspace 

Regardless of the form of the equation of state for the background, or for the perturba- 
tions, we can make the following conclusions: Any solution 

{#o, H x , H 2 , K, h , W n , W p , V n , V P7 U n , U p , 5n, dp} (75) 

to the equations governing the time-independent perturbations of a static, spherical super- 
fluid neutron star is a superposition of (i) a solution 

{H , 0, H 2 , K, 0, 0, 0, 0, 0, 0, 0, 5n, 5p} (76) 

and (ii) a solution 

{0, H x , 0, 0, ho, W n , W p , V n , V p , U n , U p , 0, 0} . (77) 

The solutions in (i) are those that satisfy the Group I equations that map one static and 
spherically symmetric star to another (nearby) static and spherically symmetric star. The 
solutions in (ii) are those that satisfy the Group II and III equations. It is not difficult to 
see that the solutions (ii) contain two sub-classes, the purely polar solutions that satisfy the 
Group II equations, i.e. 

{0, H u 0, 0, 0, W n , W p , V n , V p , 0, 0, 0, 0} (78) 

and the purely axial solutions that satisfy the Group III equations, i.e. 

{0,0,0,0,^,0,0,0,0,^,^,0,0} . (79) 

The first subclass is made of the g-modes because they are (1) purely polar and (2) the 
particle number densities (and likewise the energy density and pressure) vanish. The second 
subclass is made of the r-modes because they are (1) purely axial and (2) the particle number 
densities (and likewise the energy density and pressure) vanish. The final conclusion is 
that the zero-frequency subspace of superfluid neutron stars is spanned by two sets of g- 
modes and two sets of r-modes. This is qualitatively the same conclusion as found for 



the Newtonian equations |48| , and we assert that there will be no non-zero frequency g- 
modes in the pulsation spectrum of non-rotating superfluid neutron stars. As for the axial 
perturbations, presumably rotation will break the degeneracy and we will find two sets of 
r-modes (or more general hybrid modes [70, [71[]), as in the Newtonian superfluid case, but 
this must be verified. 
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VI. CONCLUDING REMARKS 

We have reviewed recent work to model the rotation and oscillation dynamics of Newto- 
nian and general relativistic superfluid neutron stars. We have seen that superfluidity af- 
fects both the background and the perturbation spectrum of neutron stars and both should 
therefore cause an imprint of superfluidity to be placed in the star's gravitational waves. In 
particular our local analysis of the Newtonian mode equations indicates that the superfluid 
modes should have a sensitive dependence on entrainment parameters, something that is 
supported by the quasinormal mode calculations. Given a suitably advanced detector, like 
the EURO configuration, we have seen that gravitational waves emitted during a Vela glitch, 
say, should be detectable. The key conclusion is that direct detection of gravitational waves 
from glitching pulsars can be used to greatly improve our understanding of the local state of 
matter in superfluid neutron stars. This may become more important in the next few years 
because of indications of free precession in neutron stars ||99|| , which if true means that the 
vortex creep model will have to be reconsidered. 

We have also put in place some groundwork for a future analysis of the CFS mechanism in 
superfluid neutron stars. We have done this by using the "pull-back" formalism to motivate 
fluid variations in terms of constrained Lagrangian displacements. We have also used them 
as perturbations to help map out the zero frequency subspace, and found that it is spanned 
by two sets of polar and two sets of axial perturbations. It remains to be seen if adding 
rotation will lift the degeneracy and yield two sets of polar and two sets of axial oscillation 
modes, or if the more general inertial hybrid modes result. 

Finally, I would like to elaborate a little more on why I continue to treasure my two years 
in Israel with Jacob Bekenstein. After I had completed working on the superfluid analogs of 
quantum field theory in curved spacetime effects I tried in vain to publish the work. True 
to his character, Jacob had some very kind words of advise, which were to never worry 
that effort is lost when a project does not play out exactly as expected, because his own 
experience was that one, or many, pieces of it would eventually be of direct importance for 
something else. For a young scientist, those were words of comfort and hope, and remain so 
for one that now has a little more experience. 
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Figure Captions 



Figure 1. The radial profiles of the neutron and proton background particle number densities, 
n and p, respectively, for model II. The model has been constructed such that it accords 
well with a 1.4M neutron star determined using the modern equation of state calculated 
by Akmal, Pandharipande and Ravenhall [54] . For reference, we show as horizontal lines the 



number densities at which Akmal et al suggest that i) neutron drip occurs, ii) there is an 
equal number of nuclei and neutron gas, and iii) the crust/core interface is located. It should 
be noted that the latter should not coincide with our core/envelope interface since one would 
expect there to be a region where crust nuclei are penetrated by a neutron superfluid. 

Figure 2. Comparison of one-fluid slow-rotation configurations with numerical results ob- 
tained using the LORENE code. R equ and -R po i e are the star's equatorial and polar radius 
respectively. The stellar model is a N = 1 polytrope with mass M = 1.4M and (static) 
radius R = 10 km. 

Figure 3. Typical results for the frame-dragging uj for superfluid stars in which the neutrons 
and the protons rotate at slightly different rates. 

Figure 4. The frame dragging in an extreme (rather unphysical) situation where the neutrons 
and the protons counterrotate in such a way that the net frame dragging is "backwards" at 
the surface of the star but "forwards" in the central parts. 

Figure 5. The Kepler limit Qk is shown as a function of the relative rotation rate Q n /Q p for 
our model star. The filled squares show the maximum allowed rotation rate for the protons 
(fl p ), while the open squares show the corresponding neutron spin rate (fi n ). The Kepler 
frequency simply corresponds to the largest of the two. 

Figure 6. Plots of the neutron (n) and proton (p) Kepler limits as functions of the relative 
rotation Vt n /Vt p , for a = —0.5, 0, 0.5 and e — 0, 0.4, 0.7. 

Figure 7. This figure shows the asymptotic amplitude A [n as a function of the (real) frequency 
ioM for our model I. The slowly-damped QNMs of the star show up as zeros of A{ n , i.e. deep 
minima in the figure. The first few "ordinary" and "superfluid" modes are identified in the 
figure. 
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Figure 8. This figure shows how the frequencies of the fluid pulsation modes for our model II 
vary with the entrainment parameter e. The modes shown as solid lines are such that the 
two fluids are essentially comoving in the e — ► limit, while the modes shown as dashed 
lines are countermoving. As is apparent from the data, the higher order modes exhibit 
avoided crossings as e varies. Recall that the range often taken as "physically relevant" is 
0.04 < e < 0.2. We indicate by a e and b e the particular modes for which the eigenfunctions 
are shown in Fig. ^. 

Figure 9. An illustration of the fact that the modes exchange properties during an avoided 
crossing. We consider two modes, labelled by a t and b e (cf. Figure H). The mode eigenfunc- 
tions are represented by the two Lagrangian number density variations, An and Ap (solid 
and dashed lines, respectively). For mode a the two fluids are essentially countermoving in 
the e — > limit (it is a superfluid mode), while the two fluids comove for mode b (it is similar 
to a standard p-mode). After the avoided crossing (which takes place roughly at e = 0.1) 
the two modes have exchanged properties. 
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Table Captions 



Table 1. Parameters describing our stellar models I, II and III. Model I is identical to 



model 2 of ||52|| , and has only a core with no envelope. Model II has an envelope of roughly 
1 km and could be seen as a slightly more realistic neutron star model. Model III has no 
distinct envelope but does have a more realistic mass and radius than Model I. 
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FIG. 5: 
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Tables 



TABLE I: 





model I 


model II 


model III 


a n /m n 


0.2 


0.22 


0.2 


(Jp/m n 


0.5 


1.95 


2 


Pn 


2.5 


2.01 


2.3 


P 


2.0 


2.38 


1.95 


n c (fm" 3 ) 


1.3 


1.21 


0.93 


p c (fm~ 3 ) 


0.741 


0.22 


0.095 


M/M & 


1.355 


1.37 


1.409 


R (km) 


7.92 


10.19 


10.076 


Rc (km) 




8.90 





